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ABSTRACT 

We present Gemini griz' photometry for 521 globular cluster (GC) candidates in 
a 5.5 X 5.5 arcmin field centred 3.8 arcmin to the south and 0.9 arcmin to the west of 
the centre of the giant elliptical galaxy NGC 4486. All these objects have previously 
published (C — Ti) photometry. We also present new (C — Ti) photometry for 338 
globulars, within 1.7 arcmin in galactocentric radius, which have (g — z) colours in the 
photometric system adopted by the Virgo Cluster Survey of the Advanced Camera 
for Surveys of the Hubble Space Telescope (HST). These photometric data are used 
to define a self-consistent multicolour grid (avoiding polynomial fits) and preliminar- 
ily calibrated in terms of two chemical abundance scales. The resulting multicolour 
colour-chemical abundance relations are used to test GC chemical abundance distri- 
butions. This is accomplished by modelling the ten GC colour histograms that can be 
defined in terms of the Cgriz' bands. Our results suggest that the best fit to the GC 
observed colour histograms is consistent with a genuinely bimodal chemical abundance 
distribution Nqc{Z). On the other side, each ("blue" and "red") GC subpopulation 
follows a distinct colour-colour relation. 
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1 INTRODUCTION 

Globular clusters (GCs) are tracers of early events in the star 
forming history in galaxies. However, a unique integrating 
picture of that history, beyond some tentative approaches, is 
still missing. A thorough review o f several issues in this con - 
text is presented, for example, in iBrodie fc StradeJ (|2006l l. 
Important aspects, that eventually deal with large scale 
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prop erties of galaxies (see, for example iForte. Vega fc Faifeij 
|2012| . and references therein) are both the age and chemical 
abundance distribution of these clusters. 

Even though the quality and volume of chemical 
abundance ([Z/H]) data for GCs is s teadily growing 
lAlves-Brito et al.ll20Tll : lUsher et~al]|2012l ). a key issue re- 
mains as an open subject: the connection between the GC 
abundances and their integrated colours. 

Under the common assumption of old ages, GC inte- 
grated colours should be dominated by chemical abundance 
(and in a secondary way b y age) . Evidence in t his sense can 
be found, for example, in iNorris at all ||2008D and, in the 
partic ular case of elliptical galaxies, in IChies-Santos et"al] 
("20121). 

A survey of the literature reveals numerous at- 
tempts to link colo urs and chemica l abun dance, rang- 
ing from line ar (ICeisler fc Fort3 
jHarris fc Harris! |2002|: 



1990(), q ua dratic 
Forte, Faifer fc GeisleJ l2007l : 
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dences l|Blakeslee. Cantielo fc Pendl2010l ). A recent contri- 
bution on this subject has been presented by lUsher et ahl 
l|2012l ) who adopt broken line fits. 

A clarification of the colour-abundance connection is re- 
quired, since some non linear relations (e.g. "infiected") can 
eventually lead to GC bimodal colour distributions, even 
when a unimodal chemi cal abundance di stribution is as- 
sumed (see, for example, lYoon et al.ll201ll ). Since bimodal 
colour distributions have often been identified as the re- 
sult of genuine bimodal chemical abundance distributions, 
the presence of non linearities would have important conse- 
quences on the interpretation of the GC chemical abundance 
distributions and also on their possible quantitative connec- 
tions with the diffuse stellar population in a given galaxy. 

NGC 4486 is a particularly useful galaxy in order to 
revise the chemical abundance-integrated colour issue due 
to its large GC p opulation and relative proximity to th e 
Sun, « 16.6 Mpc (|Tonrv et al.ll200ll : iBlakeslee et al.ll2009l ). 
This paper presents Gemini griz' high quality photometry 
for a selected field centred 3.9 arcmin from the centre of the 
galaxy, including 521 GC candidates, and new (C — Ti) pho- 
tometry for 338 clusters within a galactocentric radius of 1.7 
arcmin, which also have [g — z) colours obtained with the 
Advanced Camera for Surveys (ACS) of the Hubble Space 
Telescope (HST) l| Jordan et al.l l2009l ). In addition, aiming 
at extending the wavelength coverage towards the ultravi- 
olet, we combine our Gemini photometry with the C mag- 
nitudes (Wa shington system: [ Harris fc Canternalll97'ii) data 
presented bv lForte etall (|2007l . hereafter FFGO?). Ah these 
data sets can be, firstly, mutually connected to define a self 
consistent multicolour grid, and then calibrated in terms of 
different chemical abundance scales, with a final goal of ob- 
taining a simultaneous connection between metallicity and 
the colour indices grid. 

The structure of the paper is as follows. Observations 
and data handling are presented in Section[2l The relation 
between the ten colour indices defined by the Cgriz' pho- 
tometry is given in Section|3l This last section also explains 
the connection between those colour indices and others, 
like (C — Ti), {g — 2)(acs)i commonly used in extragalactic 
GCs research. Section|4]presents a preliminary calibration in 
terms of chemical abund ance using two differe nt empirical 
colour-metallicity scales: IBlakeslee et al.l ||2010|), and alter- 
natively, the scale presented by Usher et al. 1 20121 ). that are 
refreered to as BIO and U12, respectively, in what follows. 
As explained in Section^ the multicolour grid can be used 
to define different magnitude curves (or pseudo-spectral dis- 
tributions) as a function of [Z/H]. These template curves are 
used to clean the GC candidates photometric sample from 
field interlopers. The analysis of the residuals, in turn, would 
allow the detection of eventual age effects. The determina- 
tion of GCs [Z/H], via multicolor fits, is presented in Section 
[6l The results of confronting the ten observed GC colour his- 
tograms with models that adopt each chemical abundance 
calibration are described in Section!?] The final conclusions 
are given in Section[8l 




Figure 1. The GMOS field (cyan lines square; 5.5 arcmin on a 
side) discussed in this work. NGC 4486 appears up and to the left. 
North is up and East to the left. 



Table 1. Observing log of Gemini-GMOS observations. 



Date 


Filter 


Exposure 


Airmass 


Seeing 






(sec.) 




(arcsec) 


2010/01/16 


g' 


5 X 400 


1.009 


0.63 


2010/01/16-17 


r' 


5 X 300 


1.017 


0.60 


2010/01/17 


i> 


5 X 300 


1.147 


0.54 


2010/01/17 


z' 


6 X 450 


1.080 


0.47 



North) and are part of a program that also includes GMOS 
spectroscopy of a sample of selected objects that are con- 
sidered as good GCs candidates (GN-2010A-Q-21; PI: J. C. 
Forte). The field (5.5 arcmin on a side), centred 3.8 arcmin to 
the South and 0.9 arcmin to the West of NGC 4486, is shown 
in Figure[l] The observing log, including dates, filters, ex- 
posures, mean air mass and composite seeing (FWHM) is 
given in Table[T] 

Image processing was performed with the tasks of 
the package GMOS within iRAlQ. In turn, PSF magnitudes 
were obtained with the package daophot within iraf. 
These instrumental magnitudes were corrected for atmo- 
spheric extinction adopting the coefficients given in the 
Gemini web page^, and transformed to the griz' sys- 
tem using zero points derived from observations of a stan- 
dard field (PG 1323 -1-086), in a way si milar to that exten- 
sively described in iFaifer et al.l (|201ll ). The identification 
of GC candidates, using both daophot and SExtractor 



2 OBSERVATIONS AND DATA HANDLING 

The photometric observations presented in this paper were 
carried out with the Fred Gillette 8-m Telescope (Gemini 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

^ http : / / www. gemini.edu / sciops /instruments / gmos /calibration? 
q=nodc/10445 
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Table 2. Cgriz' photometry errors, in magnitudes. 



9' 






cr,/ 






20.5 


0.011 


0.010 


0.012 


0.019 


0.034 


21.5 


0.012 


0.010 


0.012 


0.020 


0.033 


22.5 


0.015 


0.011 


0.014 


0.020 


0.045 


23.5 


0.027 


0.017 


0.022 


0.035 


0.055 


24.0 


0.032 


0.021 


0.027 


0.039 


0.070 



l|Bertin fc Arnouti Il996t ) , also follows the lines described 
in that paper. The limiting magnitude of the photomet- 
ric sample is g'g « 23.7 mag, i.e., ~ 0.15 mag brighter 
than the turn-over of the GCs integrated luminosity function 
l|Villegas et al.ll2010l ). The photometric errors as a function 
of the g' magnitudes are given in Tabled 

In what follows we adopt a colour excess E(B — V) = 
0.02 m ag from the maps by ISchlegel. Finkbeiner fc DavisI 
lll998l). and the interstellar extinction relations given in 
[Jordan et al.1 (2004). The magnitudes and colours, corrected 
for interstellar reddening (denoted with the "0" subscript), 
of our GC sample are given in Table|31 

The 521 GC candidates listed in Table|3] also have 
(C — Ti) photometry obtaine d with the CTIO Blanco 4- 
m Telescope and presented in IFFGQTI. The GCs identifica- 
tion numbers are from that work. Due to severe incomplete- 
ness effects, data for GC s closer to the galaxy centre were 
not published in IFFGOTI . Among these objects we identi- 
fied 338 that wore also observed in the ACS {g — z) colour 
l| Jordan et al.i ,2009). The photometric values for these GC 
candidates are listed in Table|3]where the R.A. and Dec val- 
ues come from this last work. 

The distribution on the sky for objects within our Gem- 
ini field are depicted in Figure[2l while their go vs. {g — z)'q 
and Tiq vs. (C — Ti)o colour magnitude diagrams are shown 
in FigureO 



3 MULTICOLOUR RELATIONS 

The adoption of a polynomial fit to represent the depen- 
dence of integrated GC chemical abundance \Z/H] with 
some particular colour index (or vice versa), is a usual ap- 
proach in t he literature. For example, a linear fit was at- 
tempted by iGeisler fc Forti (|l990l ). while two second or- 
der polynomials, minimizing t he errors in ab u ndanc e or 
in colour, were presented by Harris fc If arri^ 1*2 0021). In 
turn, both iMovano Loyola et alT ( 2010l ) aiid Blakcs lee et al.l 
l|2010t ) adopted the robust routine, in their respective sec- 
ond and quartic order approximations. That routine seeks 
minimizing the "orthogonal" errors (defined in terms of a 
residual that combines those of the two variables). This ap- 
proach was also adopted by Blakcslee et al. (2012) when de- 
riving the relation between the {g — I) and (/ — H) colours 
of GC candidates in the central field of NGC 1399. 

In this paper we attempt a different approach aiming 
at connecting simultaneously the ten different colour in- 
dices that can be defined through the Cgriz' photometry 
and avoiding polynomial fits. The resulting self-consistent 
colour-grid is then tentatively calibrated in terms of two 
chemical abundance scales, as explained in the following sec- 
tion. 
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Figure 2. Distribution on the sky of 521 GC candidates with 
Cgriz' photometry Hsted in Table[3] North is up, East to the left. 
The centre of NGC 4486 is at a; = 0.0, y = 0.0. 
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Figure 3. vs. {g — z)g and Tig vs. (C — Ti)o colour magnitude 
diagrams (corrected for interstellar reddening) for the GC can- 
didates listed in Table[3] Note incompleteness effects in the first 
panel for clusters bluer than [g — z)'q = 1.25. 



As a first step we generated nine colour-colour planes, 
all including the {g — z)'q colour, the index with the longest 
wavelength base (in terms of the Gemini photometry). We 
preferred that index instead of (C — z')o given the larger er- 
rors inherent to the C magnitudes (see Table[2]). Each pho- 
tometric value was then convolved with a Gaussian kernel 
to generate bi-dimensional images that were afterwards used 
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Table 3. Multicolour photometry for GC cand idates in the NGC4486 field. Magnitudes and colours are corrected for interstellar extinction 
and reddening. Identification numbers are from lpbrte et al" ] {2007). The full version of the table will be available in the electronic edition of 
the Journal. 



ID 


9'o (C-9')o (C-r')o (C-i')o (C-^')o 


{9-r)'o 


{9 - 


(9 - ^)(, 






- ^);, 


{C-Ti)o 


23314 
23796 


22.38 0.62 1.26 1.50 1.61 
21.71 0.38 0.93 1.14 1.19 


0.64 
0.56 


0.87 
0.76 


0.99 
0.81 


0.23 
0.21 


0.35 
0.25 


0.12 

0.05 


21.62 1.39 
20.95 1.14 



Table 4. GC candidates with ACS and Washington Photometry within a galactocentric radius of 100 arcsecs in NGC 4486. The full 
version of the table will be available in the electronic edition of the Journal. 



a 


<5 




90 


(C-Ti)o 


(9 - 2^)0(ACS) 


(degrees) 


(degrees) 


(mag) 


(mag) 


(mag) 


(mag) 


187.7037354 


12.3645134 


22.686 


23.429 


1.140 


0.900 


187.7002411 


12.3649130 


22.876 


23.701 


0.940 


0.789 




Figure 4. Smoothed colour-colour relations. From top (g — r)^, 
(1 — 2)q, and {i — z)^ (shifted arbitrarily in ordinates) vs. {g — z)g. 
The frame is 1.1 mag on a side. 



to obtain modal colour values. The adopted kernel size (0.05 
mag) matches the overall error of the colours in our sample. 

These images, generated and processed with iraf 
(through the routines irafil and gaussian), were used to 
determine the modal values linking each index at a given 
((? — z)'q colour (from 0.80 to 1.50 with 0.1 mag intervals, 
and adding two points at the bluest and reddest ends: 0.75 
and 1.55 mag, respectively). Outside these limits the de- 
termination of the modes becomes uncertain. An example, 
corresponding to the ((/ — r)o, (j" — i)o s^nd (i — z)o vs. (g — 2)0 
relations is depicted in Figure[4l 

We also performed a number of simulations start- 
ing with different types of colour-colour relations (linear, 
quadratic, fourth order) and adding Gaussian errors (simi- 
lar to those in Table 2). The smoothing procedure was able 
to recover the input relation with a typical uncertainty rang- 
ing from ± 0.012 mag to ± 0.02 mag. These simulations also 
show that changing the smoothing kernel from 0.03 to 0.08 
mag has no detectable impact on the derived modal colours. 



In fact, each modal colour value can be determined in 
four different ways. First, directly from the plane determined 
by a given colour vs. (g — 2)0 and then by properly combining 
all the other indices including the filter bands that define 
that colour (e.g., the {C—g')o vs. ((7—2)0 relation can also be 
determined from the [(C — r')o — (<? — r)o], [(C — i')o — (3 — 1)0], 
and [(C — z')o — ((? — z)o] modal colour differences determined 
at a given [g — z)'q). These four indices will differ as a result 
of the distinct roles played by the photometric errors on each 
of the two colour planes. 

A check for consistency indicates that, for a given mean 
colour, the scatter of each four modal colours has a typical 
rms of « 0.012 mag (and « 0.018 mag in the worst case) 
over the whole grid. The link to the (C — ri)o colour was 
obtained adopting the same method described above, i.e., by 
generating two smoothed colour planes and looking for the 
modal colour values. In turn, (C — ri)o was also connected 
to (g — z)o^Apg-, thr ough the GC candidates in common with 



[Jordan et 411(20091 ) . listed in Table g] 
From the colour grid we derive: 

(C - ri)o = 1.25 (±0.03) {g ~ 2)o(acs) + 0-08 (±0.04). (1) 

This equation corresponds to the formal bisector fit. We 
stress, however, that there is a deviation at the blue extreme 
that, as discussed later, seems to be a common feature in all 
the colour-colour relations involving the C filter, possibly as 
a consequence of the presence of the features connected with 
the 4000 A break. From the same grid, we obtain: 



(5 - = {g- ^)o(Acs) - 0.07 (±0.02). 



(2) 



The relations between (C — ri)o and [g — z)o with 
((? — -2)0 (ACS) gi"^6n in Table[5]are displayed in Figure[5] 

The colour-colour relations defined by Table[S] display 
different degrees of non linearity. In some cases, these effects 
can be noticed near the blue and red ends of the colour- 
colour relat ions. Similar kind of no n linearities were already 
noticed by i Blakeslee et al.l (|2012| ) in their analysis of the 
(optical) [g — z)acb vs. the (infrared) (J — H) index. As 
examples, Figure[B] and Figure[7] show the behaviour of all 
the colour indices defined in our photometry as a function 
of (C — g')o and (g — z)o, respectively. 
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Table 5. Multicolour-Chemical abundance relations from GCs in the field of NGC4486 



{C -g')o 


(C - r')o 


(C-i')o 


(C ~ z')o 


iq-r)',. 




{g - zYg 


(r - iYn 


(r ~ zYn 


(i - zYn 


(C - Ti)o 


{9 ~ z)o 
ACS 


[Z/H] 
BIO 


[Z/H] 
U12 


0.29 


0.82 


0.99 


1.04 


0.52 


0.70 


0.75 


0.18 


0.23 


0.05 


1.05 


0.83 


-2.16 


-1.87 


0.38 


0.93 


1.11 


1.18 


0.54 


0.73 


0.80 


0.19 


0.26 


0.07 


1.15 


0.86 


-1.88 


-1.64 


0.51 


1.10 


1.32 


1.41 


0.59 


0.80 


0.90 


0.21 


0.30 


0.09 


1.27 


0.94 


-1.35 


-1.24 


0.60 


1.24 


1.49 


1.60 


0.63 


0.88 


1.00 


0.25 


0.37 


0.12 


1.42 


1.06 


-0.84 


-0.96 


0.69 


1.36 


1.64 


1.79 


0.68 


0.96 


1.10 


0.28 


0.43 


0.14 


1.53 


1.14 


-0.64 


-0.68 


0.76 


1.48 


1.79 


1.96 


0.72 


1.03 


1.20 


0.31 


0.48 


0.17 


1.67 


1.26 


-0.44 


-0.42 


0.86 


1.61 


1.95 


2.16 


0.74 


1.09 


1.30 


0.35 


0.56 


0.21 


1.81 


1.38 


-0.23 


-0.22 


0.96 


1.74 


2.11 


2.36 


0.79 


1.16 


1.40 


0.37 


0.61 


0.24 


1.91 


1.48 


0.09 


0.02 


1.04 


1.86 


2.25 


2.54 


0.82 


1.21 


1.50 


0.39 


0.68 


0.29 


2.04 


1.59 


0.70: 


0.19 


1.08 


1.92 


2.33 


2.63 


0.85 


1.25 


1.55 


0.40 


0.70 


0.30 


2.11 


1.64 




0.33 



I 

00 





(g-z)„ 



(ACS) 



(c-g')„ 



Figure 5. From top to bottom: {C — Ti)o and (g — z)^ vs. 
{g — z)o(ACS) relations. Filled dots are the modal colours listed 
in Table0 The straight lines correspond to Equation[l]and Equa- 
tionlll 



Figure 6. Colour indices vs. (C — g')o- Prom top to bottom, 
in decreasing order: [C — z')o, (C — i')Oi (C* — r')o (open dots); 
(9 - 2);,, {g - «);,, (g - tYq (filled dots); (r - 2)^, (r - i)^, (open 
squares); (i — z)'q (filled squares). The straight line at upper left 
has a unity slope. 



4 COLOUR-CHEMICAL ABUNDANCE 
RELATIONS 

Despite significant effort, the shape of the relation between 
GC integrated colours and chemical abundance still remains 
a subject of debate. The main problems behind a proper 
determination of their relation are both the uncertainties of 
synthetic models and the still large errors associated with 
the GC che mical abundances deriv ed, for example, via Lick 
indices fsee lBrodie Sz Huchralll99ol '). 

As a preliminary approach, and in order to assess how 
the inferred abundances depend on the adopted calibration, 
we use two different empirical calibrations. On the one hand, 
we adopt the "inflected" (g — z) acs vs. [Fe/H] relation pre- 
sented by BIO and, on the other, the "broken line" {g — i) 
vs [Z/H] relation given by U12. The first one relies mostly 
on Galactic GCs and includes some high [Z/H] globulars 
from NGC4486 and NGC4472. In turn, the second relation 



stands on a compilation of the literature and spectroscopic 
observations of the Calcium triplet lines. 

We remark that the upper abundance value given in 
Tables, corresponding to the BIO calibration, is just a for- 
mal extrapolation and, in what follows (e.g. model fits), we 
adopt an upper cut-off at [Z/H] = 0.40. The same comment 
holds for the U12 calibration that only reaches solar abun- 
dances and is extrapolated up to [Z/H] = 0.33. With this 
caveat in mind, we note that only a small fraction of the GC 
candidates in our sample seem to have abundances higher 
than solar. 



5 GLOBULAR CLUSTER CANDIDATES 

In Figure |8] and |9] we display the ten colour histograms de- 
fined in terms of the Cgriz' photometry presented in this 
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Figure 7. Colour indices vs. {g — 2)g. From top to bottom, in 
decreasing order: (C — z')o, {C — i)g, (C — r-')o, (C — g')o (open 
dots); (g - i)^, (g - r)[, (filled dots); (r - z)'^, (r - i)'^ (open 
squares); {i — z)q (filled squares). The straight line at upper left 
has a unity slope. 



Figure 9. Colour histograms (from top to bottom, left panel): 
{g — r)'o, (g — i)'o, (g — z)' o and (right panel, from top to bottom) 
(r — i)'o, (r — z)'o and (i — 2:)'o for all the objects listed in Table 
3. The histograms have beeen shifted arbitrarily in ordinates. 
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Figure 8. Colour histograms (from top to bottom): (C — g')o, 
{C - r')o, (C - i')o, {C - z')o for all the objects listed in Table 
3. The histograms have been shifted arbitrarily in ordinates. 
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Figure 10. Magnitude curves defined by the modal colours listed 
in TableOfrom (g — z)q = 0.80 to 1.50 (in 0.1 mag steps; contin- 
uous lines). The dashed upper and lower lines correspond to the 
extreme colours (g — z)^ = 0.75 and 1.55, respectively. 



paper. These colour distributions correspond to all the ob- 
jects listed in Table[3] This sample includes a fraction of field 
interlopers. In order to decrease the effect of these objects in 
our following analysis, we use the template values listed in 
Table [5] in an attempt to identify the genuine GCs. We stress 
that the rejected objects, however, have a negligible effect 



on the definition of the modal colour-colour relations. That 
table defines, at each [Z///], a template magnitude curve (or 
pseudo spectral distribution) as shown in Figure fTOl In this 
diagram we arbitrarily adopt a reference magnitude Zq = 0.0 
mag. These curves assume that GCs are coeval on the basis 
of the arguments given before. 
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Figure 11. Root mean square values corresponding to the mag- 
nitude curve fits as a function of tlie i'^ magnitudes for tlie 521 
objects included in Table 3 and adopting the U12 chemical abun- 
dance scale. Objects belowr 0.05 mag (dotted line) are considered 
as GC candidates. 



Figure 12. Upper panel: Magnitude residuals {C, g' ,r' ,i' , z' 
bands) for the whole sample (521 objects). Dashed horizontal 
lines indicate the adopted rms = ±0.05 mag limits that define 
the "clean" GC sample. Lower panel: Expected residuals from 
Monte Carlo modelling for an age spread of ±2Gyr (see text). 



The magnitude curves were used to fit the {Cgriz')o 
magnitudes of all the objects listed in TableO i.e., we vary 
the [Z/H] value by interpolating in Table[5] and adopting 
0.01 steps in [Z/H] until a given template curve minimizes 
the sum of the square residuals at each filter band. This pro- 
cedure was performed using each of the chemical abundance 
scales. 

As the results are almost identical in terms of accepted 
and rejected GC candidates, and for illustrative purposes, 
in this Section we only show the diagrams corresponding to 
the adoption of the U12 calibration. 

Figure fm shows the overall composite rms (defined by 
combining all the magnitude residuals) as a function of the 
if) magnitude. We adopted this magnitude, in particular, be- 
cause in this way all the photometric errors in that diagram 
are decoupled. The dotted line indicates a rms = 0.05 mag 
that we adopt as a "reasonable" maximum acceptable value 
to consider a given object as a genuine GC candidate. With 
this criterion, the BIO scale accepts 463 objects as GC can- 
didates and rejects 58, while the U12 scale leads to 472 and 
49, respectively. The difference in the number of rejected ob- 
jects arises as a consequence of the upper colour limit in the 
BIO calibration (see Table 5). This limit produces rms val- 
ues larger than the adopted rejection limit for nine objects 
with colours redder than [C — g')o ~ 0.95, (r — 2)0 — 0.65, 
and (C - z')o = 2.5. 

The upper panel in Figure[T2]shows the residuals yielded 
by this procedure at each filter (characterized by their effec- 
tive wavelength). These residuals do not show any system- 
atic trend. 

As a test of the effects due to an eventual range in the 
GC ages, we generated a Monte Carlo model that includes a 
spread of ±2 Gyr on top of the colour-abundance relations 



given in TableO The differential magnit ude variat i ons a s 
a function of age were obtained using the iMarastonI (120051 ) 
models and adopting the 12 Gyr models as the reference age. 

The output colours, including age variations but not 
photometric errors, were then fit with the GC templates, as 
in the case of the observed clusters, yielding the residuals 
displayed in the lower panel of Figure fT^ The behaviour of 
these residuals indicates that age effects have no important 
impact on the fits as they are considerably smaller than 
the photometric errors. On the other hand, the same Monte 
Carlo models show that, with our photometric errors, the 
input [Z/H] values can be recovered through the magnitude 
curve fit with an rms « 0.15 dex. 

The number of rejected objects is consistent with previ- 
ous results based on CMOS spectroscopy. In this last case, 
the photometric criteria adopted to define a GC candidate 
typically yield a contamination level of 10 perce nt once these 
cand idates are observed spectroscopically (see iFaifer et al.l 
2011i). 

The residuals of the (C — g')o and (r — «)o colours as 
a function of the i'g magnitude are shown in Figure[T31 GC 
candidates fainter than ig ~ 20.5 mag do not show obvi- 
ous systematic behaviours. However, small anti-correlated 
trends are observed for the brightest objects. 

We have been unable to reproduce the residual trends 
seen for the brightest objects just by changing the cluster 
ages in our models. Neither the photometric errors, nor those 
connected with the inferred shape of the multicolour rela- 
tions, seem adequate explanations for this effect. 

The clean GC sample, and the rejected objects, are 
shown in Figures [T4l and [151 where we transformed {g — zYq 
to {g — •z)o(ACS) through Equation^ Both diagrams show 
non linear behaviour with detectable changes in their slopes 
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Figure 13. Fit residuals in the (C — g) and (r — z)q colours as a 
function of «q magnitude for all the GC candidates adopting the 
U12 calibration. 

at (C — g')o ~ 0.5 and (C — z')o « 1.4, respectively. In 
particular, Fisure |1 5 | can be compared with figure 16 in 
IChies-Santos et all j201ll ). which displays (<? — z)o(acs) ^ 
a function of the optical-infrared [g — K) index for GCs in 
NGC 4486 and NGC 4649. Our own reading of that diagram 
indicates that GCs in these galaxies show changes in the 
colour-colour slope at {g — K) fa 2.90 and {g — K) ~ 3.5. 
In turn, GCs bluer than the first of these colours span a 
{g — z)acs range that covers from « 0.8 to « 1.0 , i.e., the 
region where we detect a change in the colour-colour slope 
in Figure lTS] 



6 INFERRING THE GC [Z/H]S THROUGH 
THE TEMPLATE PSEUDO-CONTINUUMS 

In this section we make an attempt to recover the GCs [Z/H] 
distribution using their integrated colours. A simi l ar ap - 
proach is presented, for example, bv lBlakeslee et al.1 (|2012l ') 
(see their figure 11, right panel) for the case of the NGC 1399 
clusters. Instead of using a single integrated colour, in this 
work we obtain the [Z/H] values for the GC candidates from 
the "template curve" fitting described in the previous sec- 
tion. We infer two chemical abundance distributions by al- 
ternatively adopting the BIO or U12 calibrations. 

As discussed in several papers, the GC colour bimodal- 
ity in NGC 4486 is better define d for G Cs fainter than 
Tio « 21.0 mag (see, for example. [fFGOTI or iHarrid [2OO9I ') 
and, accordingly, we split our sample in two groups: GCs 
candidates with Ti from 19.0 to 21.0 mag and Ti from 21.0 
to 23.0 mag (i.e. ~0.2 mag brighter than the turn-over of 
the integrated GC luminosity function). 

Figure [T6l and FigurelTTi exhibit the resulting [Z/H] dis- 
tributions adopting the BIO and U12 calibrations respec- 
tively and for the two GC groups. In both diagrams the 
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Figure 14. {r — z)'q vs. {C — g')o colours showing accepted (filled 
dots) and rejected (open dots) CG candidates adopting the U12 
calibration and after applying the rms > 0.05 mag cleaning crite- 
rion. The straight line at upper left has a unity slope. The lower 
curve corresponds to colours given in TableOshifted by —0.5 mag 
in (r - 2)^. 




(C-z')„ 

Figure 15. [g — z)o(aqs) ~ ^')o colours for the accepted 

(filled dots) and rejected (open dots) GC candidates adopting the 
U12 calibration. The straight line at upper left has a unity slope. 
The lower curve corresponds to colours given in Table[5] shifted 
by -0.5 mag in {g - z)o(acS)- 
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brightest clusters show distinct behaviour compared with 
the fainter counterparts. For the brighter GCs the abun- 
dance distributions seem broad and unimodal. 

In turn, independent of the cahbration, GCs fainter 
than Ti — 21 mag show a clear bimodality although the 
histograms exhibit differences in their shapes and in the po- 
sition of the low and high abundance peaks (that do depend 
on the adopted calibration). 

A simple Gaussian analysis (using the RMIX routine) 
indicates that, for these clusters, and for both calibrations, 
two Gaussian fits are strongly preferred over a single Gaus- 
sian fit, leading to values about 3 to 5 times smaller. 
For the case of the BIO calibration {x^ = 1.87) we obtain 
mean [Z/H] values of —1.66 and —0.37 with dispersions of 
0.32 and 0.39 dex for the blue (55 percent of the population) 
and red GCs, respectively. Alternatively, adopting the U12 
calibration (x^ = 0.95), these parameters become [Z/H] of 
— 1.42 and —0.28, with dispersions of 0.29 and 0.30 dex for 
the blue GCs (61 percent of the population) and red GC 
populations, respectively. 

As these Gaussians have comparable dispersions, we 
also attempted a homoscedastic KMM test that, in both 
cases, indicates that the probability for a single Gaussian fit 
to represent the [Z/H] distributions is practically null. 

The situation of GCs brighter than Ti — 21.0 mag is 
not so clear. Both the BIO and U12 calibrations lead to 
broad unimodal distributions that show a different degree 
of skewness. 

The results shown in Figure[T6] and Figure[T7| are differ- 
ent, but not ne c essari ly in conflict, with those presented by 
iBlakeslee et"all (l2012h for the central region of NGC 1399 
where red GCs are clearly the dominant subpopulation. 
These authors find a broad unimodal distribution with a sin- 
gle peak (—0.3 dex) and a tail extending towards lower chem- 
ical abundances (see their figure 11, right panel). In turn, 
our GCs field is located at a larger galactocentric distance, 
where the relative number of blue globulars is considerably 
larger making more evident the presence of a low chemical 
abundance peak. 



7 MODELLING THE GLOBULAR CLUSTER 
COLOUR HISTOGRAMS 

Two key issues in modelling the GC colour histograms are 
the adopted integrated colours vs. chemical abundance cali- 
bration and the assumed chemical abundance distribution of 
the clusters. In general, most papers have relied on a single 
colour-abundance calibration and also on a single GC colour 
histogram fit. In this work, we attempt a simultaneous fit to 
the ten colour histograms that can be defined in terms of the 
Cgriz' magnitudes,^adopting the (inverse) quality fit indica- 
tor x^, given in ICote. Marzke fc WestI ([1991). Finally, we 
identify the best fit parameters (that define the GC chemi- 
cal abundance distribution) as those that yield a minimum 
value for the sum of the individual indices of the ten 
colour histograms (adopting the same colour bin for all the 
histograms: 0.10 mag). 

Monte Carlo model histogra ms were obta i ned fo llowing 
the same procedure explained in iForte et al.l (|2012t ). First, 
we generate "seed" globulars with chemical abundances Z 




-2 -1 
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Figure 16. Inferred chemieal abundance distributions adopting 
the BIO calibration for GC candidates brighter than Ti = 21.0 
mag (shifted upwards by 130 units) and for GC candidates with 
Ti from 21.0 to 23 mag. Open dots correspond to the formal 
counting uncertainties for each bin. 
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Figure 17. Inferred chemical abundance distributions adopting 
the U12 calibration for GC candidates brighter than Ti = 21.0 
mag (shifted upwards by 130 units) and for GC candidates with 
Ti from 21.0 to 23 mag. Open dots correspond to the formal 
counting uncertainties for each bin. 
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(within the range Zi to ^max) whose numbers are controlled 
by a given distribution function 7Vgc(-^)- 

After trying different simple distributions, we con- 
cluded, as in FFG07, that a double exponential dependence: 



iVGc(^) « e^p{~Z/Zs 



(3) 



(where Zs is the scale length corresponding to the blue or 
red GC subpopulation) is the simplest function that allows 
a fit to the colour histograms based on a minimum number 
of free parameters. Formally, this approach requires seven 
parameters: the ratio of blue to red clusters, and for each 
GC subpopulation the Zi and Zmax values as well as the 
chemical scale lenght Zs. In fact, the free parameters were 
reduced to five, as the minimum chemical abundance of the 
blue GCs subpopulation, as well as the maximum chemical 
abundance of the red GCs, were set as the lowest and upper 
values in the BIO and U12 calibrations. 

The integrated colours for each GC were obtained by 
linear interpolation, using the logarithmic abundance as ar- 
gument, in Table[51 

For each synthetic cluster we also generate an appar- 
ent magnitude qq, adopting a Gaussian integrated lumi- 
nosity function and according to the parameters given by 
IVillegas et al](|2010l ). These magnitudes were used as input 
in Table[2]in order to model (also Gaussian) observing errors 
that were added to each colour. Given the relatively short 
range in apparent magnitude, we do not include an explicit 
dependence of chemical abundance with brightness for the 
blue globulars (i.e., the "blue tilt" effect; see, for example, 
[Harr is 2009). 

The parameters that define the chemical abundance dis- 
tributions and provide the best overall fit to the ten colour 
histograms in each case, are listed in Table[Sl and the corre- 
sponding individual and cumulative indices are given in 
TableQand Table[Hl 

Even though both models lead to N{GC) values within 
« 1.5 times the formal counting errors of each histogram bin, 
the U12 calibration yields better fits in terms of the cumu- 
lative quality index. We remark that this statement is valid 
only if the bi-exponential parametrization of the chemical 
abundance is accepted. 

Each of the parameters listed in Table[n]has an associ- 
ated uncertainty which we define as the parameter variation 
that leads to a decrease of the fit quality, indicated by an 
increase of ten percent above the mimimum total value 
in the five free parameters space. 

Following this, we get typical uncertainties of ± 10 GCs 
for each subpopulation; for the blue GCs: ± 0.01 Zq in Zsb 
and ± 0.05 Zq in Zmax; and for the red GCs: ± 0.06 Zq in 
ZsR and ± 0.02 Zq in Zi. 

The histograms corresponding to (C — z')o colour are 
depicted in Figure[T8]and Figure fTOl We note that this index 
is « 1.6 and ~ 2.0 times more sensitive to metallicity than 
(C — Ti) and {g — z'), respectively. This allows the adop- 
tion of a larger colour bin (0.2 mag) thereby reducing the 
sampling noise. Both histograms exhibit a colour "valley" 
at (C — z')o ~ 1.7, i.e., 0.3 mag. redder than the colour 
where a change in the colour-colour relation is detectable 
(see Sections). The models show that both GC populations 
overlap in the color range (C — z')o ~ 1.30 — 1.70. In this 
range, some 20 to 25 percent of the total number of clus- 
ters belong to the "blue tail" of the red GC population. We 
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Figure 18. GC (C — z')o colour distributions adopting the BIO 
calibration (solid lines) for GC candidates brighter than Ti = 21.0 
mag (shifted upwards by 130 units) and for GC candidates with 
Ti from 21.0 to 23 mag. Open dots correspond to the formal 
counting uncertainties for each bin. Dashed lines correspond to 
the best fit models given in Table 6. 



Table 6. Parameters of the chemical abundance distributions 
giving the best overall fit to the colour histograms for Globular 
Clusters fainter than Ti = 21. 
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note that the (C — z')o colour of the valley corresponds to 
(q — 2)n,Ar^o, ~ 1.1 to 1.2, i.e., coincident with the colour 
region where IChies-Santos et al.l (12011! ) claim the presence 
of a "wavy feature" in the {g — 2;)acs vs. {g — K) relation. 

The GCs [Z/H] distributions inferred from the photo- 
metric observations can be compared with those arising in 
the models that deliver the best fits to the ten observed 
colour histograms and whose parameters are listed in Ta- 
ble El 

Figures [20l and [211 displav such a comparison. For clus- 
ters fainter than Ti = 21.0 mag the models based on the 
BIO calibration yield mean [Z/H] of -1.69 and -0.47 while 
the adoption of the U12 calibration gives somewhat larger 
values, —1.46 and —0.36. These mean abundances are close 
to, but different from, the values arising from the simple 
Gaussian analysis presented before. 

The fits to the brighter clusters lead to an am- 
biguous situation given the small number of GCs 
and the absence of definite peaks that prevents us 
from obtaining robust results. 
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Table 7. GC colour histograms fit adopting the BIO calibration. First line, GCs brighter than Ti = 21.0 mag; second line: GCs fainter 
than Ti = 21.0 mag. 



{C-g)'o {C-rW 


(C-i)o' 


{C-z)o' 


(9 - rYo 


{9-i)o' 


(9 - ^Yo 


(r - iYo 


(r - ^Yo 


- zYo 


Cumulative x'^ 


0.85 0.64 


0.90 


0.74 


0.61 


0.92 


0.41 


0.08 


1.41 


1.42 


7.98 


0.56 1.18 


1.07 


1.13 


0.84 


1.29 


1.62 


0.41 


1.76 


1.60 


11.46 



Table 8. GC colour histograms fit adopting the U12 calibration. First line, GCs brighter than Ti = 21.0 mag; second line: GCs fainter 
than Ti =21.0 mag. 
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Figure 19. GC (C — z')o colour distributions adopting the U12 
calibration (solid lines) for GC candidates brighter than Ti = 21.0 
mag (shifted upwards by 130 units) and for GC candidates with 
Ti from 21.0 to 23 mag. Open dots correspond to the formal 
counting uncertainties for each bin. Dashed lines correspond to 
the best fit models given in Table 6. 



Figure 20. Comparison between the inferred [Z/H] GC distri- 
butions (solid line histogram) for GCs fainter than Ti = 21 mag, 
and that derived from the colour histogram fits adopting the BIO 
calibration (dashed lines). Open dots are the formal counting un- 
certainties in each bin. 



8 CONCLUSIONS 

This paper presents a self consistent multicolour Cgriz' grid 
including 100 points, each with a typical uncertainty of 
±0.012 mag, and based on GC candidates in NGC 4486. 
This grid, once calibrated in terms of two different colour- 
metallicity relations, has been used to infer the GCs chem- 
ical abundances from Cgriz' photometric data and to per- 
form a comparison with simple models. 
The main results are: 

1) The multicolour Cgriz' relations show different degrees 
of non linearities. This is more evident in those colours in- 
volving the C filter. Non linearities o f this kind had been 
previously reported bv lBlakeslee et all ()2012h in their study 
of the central regions of NGC 1399. 



2) The inferred [Z/H] distributions for GCs fainter than 
Ti = 21.0 mag are bimodal either adopting the "inflected" 
BIO or the "broken line" U12 colour-abundance calibrations. 

3) The model fits based on a double exponential depen- 
dence of the number of clusters with chemical abundance Z 
both for the blue and red GCs provide a good representation 
of the GC integrated colour histograms and of their inferred 
chemical abundance distributions. 

4) For the brightest clusters the abundance distributions 
appear broad and skewed and we do not reach a definite con- 
clusion regarding the presence of bi-modality. These objects 
leave systematic colour residuals from the template magni- 
tude curves that cannot be easily accounted for. In a spec- 
ulative way, these residuals might indicate the presence of 
multi-stellar populations similar to those found in systems 
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Figure 21. Comparison between the inferred [Z/H] GC distri- 
butions (solid line histogram) for GCs fainter than T\ = 21 mag, 
and that derived from the colour histogram fits adopting the U12 
calibration (dashed lines). Open dots are the formal counting un- 
certainties in each bin. 

with comparable absolute magnitudes (e.g. Ide Boer et al.l 
|2012| ). 

5) Adequate two-colour diagrams, such as (C — g')o vs. 
(r — z)'q or (C — z')o vs. {g — z)'q, show changes in the colour- 
colour slopes. These changes are detectable, for example, at 
{C ~ z')o — 1.4. On the other hand, the bi-exponential mod- 
elling, adopting the BIO or U12 colour-chemical abundance 
relations, show that 88 and 82 percent, respectively, of the 
"blue" GCs are bluer than that colour. This indicates that 
the "blue" and "red" GC sub-populations exhibit distinct 
colour-colour relations with a transition zone possibly be- 
tween (C - z')o = 1.40 and (C - z')o = 1.60 (where both 
sub-populations overlap to some degree). 

The results presented in this paper suggest that the origin 
of the GC colour bimodality has its roots in a real bimodality 
of their chemical abundance distributions, and are consistent 
with the spectroscopic analysis of GCs in roughly half of 
the sample of eigh t galaxies studied by U12 and also by 
iBrodie et. al (|2012l ) in NGC3115. 
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